function diff = v_fun( u, lambda, param)

    % Vacancies 
    diff= (lambda./param.A).^(1./(1-param.epsilon)).*(param.s*(1-u) + u);

end

